Magnetotunneling spectroscopy as a probe for pairing symmetry determination in 

quasi-2D anisotropic superconductors 
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As a probe to determine the pairing symmetry of quasi-two-dimensional anisotropic superconduc- 
tors, we propose tunneling spectroscopy in the presence of magnetic field, where the magnetic field 
is parallel to the two dimensional planes and rotated. As a case study, we apply this idea to the 
models of high-Tc cuprates and organic superconductors k-(ET)2X. The surface density of states 
at the Fermi energy exhibits a characteristic oscillation upon rotating the direction of the magnetic 
field due to the Doppler shift of the energy of quasiparticles. The surface density of states has a 
minimum when the applied magnetic field is parallel to the node direction of the pair potential in- 
dependent of the detailed shape of the Fermi surface. The amplitude of the oscillation is sensitively 
affected by the shape of Fermi surface. 



Nowadays, there are various superconductors in which 
the pairing mechanism seems to be unconventional. In 
order to clarify the pairing mechanism of these super- 
conductors, it is crucial to identify the pairing sym- 
metry, which is characterized by the sign change and 
the presence of nodes in the pair potential. In or- 
der to determine the pairing synMHstry, several phase- 
sensitive probes have been usedula^u, among which is 
the tunneling spectroscopy, which enables us to detect 
the sign change_ijji the pair potential as well as its 
nodal structureQEra. Namely, a reliable evidence of the 
sign change of the pair potential is aa observation of 
a zero-bias conductance peak (ZBCP)l3, which is orig- 
inated from Andreev bound., states (ABS) formed at 
the surfaces or interfacesBUHatl. The existence of ZBCP 

has been a r-tn a 11 y Obp^iV^iy^ ^'.Y M '^'{^'V." l| ^ ''rirrimrmiSnn a 1 

superconductorsBMaMeyyyy Bma . 

In principle, it is possible to determine the pairing sym- 
metry through the tunneling spectroscopy via ABS if one 
can make well oriented surfaces or interfaces for arbi- 
trary orientations. However, in quasi- two dimensional 
(2D) superconductors, it is by no means easy to make 
a well oriented surface/interface in the direction perpen- 
dicular to the 2D planes. For example, for the high- 
Tc cuprates, it is predicted theoretically that ZBCP is 
observed most prominently for (110) surface/interface, 
and is not observed for (100) oriented surface/interface.cl 
However, in the actual experiments, it is not easy to 
prepare well contmllfjd- surface/interface, and only few 
exceptional caseal3llj'E£l have succeeded in discriminat- 
ing the difference between the tunneling conductance of 
(100) and (110) oriented junctions. Moreover, it has been 
clarified that atomic-sqale-muahness has influence on the 
tunneling spectroscopycHijEj'Eil. Thus, in order to clearly 
determine the position of the nodes in the pair potential. 
I.e., the symmetry of the pair potential, through tunnel- 
ing spectroscopy via ABS, a preparation of well oriented 
junctions is necessary. 

The situation can even be more severe for quasi- 



2D organic superconductors such as k-(BEDT-TTF)2^ 
saltsEll [abbreviated k-(ET)2X-, hereafter], {X = 



)2fn 

Cu(NCS)2, Cu[N(CN)2]Br, or l^ifA, which is another can- 
didate for anisotropic pairing superconductors. It is even 
more difficult to make well oriented junctions compared 
to the high-Tc cuprates. 

For quasi-2D superconductors, it is much more promis- 
ing to make films well oriented in the direction parallel to 
the planes, but the problem is that in that case, ZBCP 
is not obtained in the tunneling spectroscopy due to the 
absence of the interference of the sign change of the pair 
potential felt by the quasiparticle. Thus, even if there 
are nodes in the pair potential, only a V-shaped tun- 
neling conductance is expected, and the direction of the 
nodes cannot be determined. In order to overcome this 
difficulty, we propose in the present study tunneling spec- 
troscopy for surfaces parallel to the planes in the presence 
of a magnetic field applied parallel to the planes. 

In the presence of a magnetic field, it is known that 
the energy spectrum j^it, the quasiparticle is influenced 
by the Doppler shiftBES. Due to this effect, speciflc 
heat and thermal conductivity show characteristic os- 
cillation by rotating the direction of the magnetic field 
within the 2Djilane.Ee|lecting the nodal structure of the 
pair potentia£3t3'Ej'E3. Recently, based on this idea, 
pairing symmetry of several unconventional supercon- 
ductors have becH-j disejissed using thermal conductiv- 
ity measurementa2a'E3'E3. Although thermal conductiv- 
ity and speciflc heat measurements are powerful ways 
to determine the pairing symmetry, these quantities are 
strongly influenced by phonons, so that the actual anal- 
ysis of the experimental data, namely the subtraction 
of the phonon contributions, can be subtle. Thus, it is 
desirable to propose a complementary method to deter- 
mine the positions of nodes. In this sense, tunneling spec- 
troscopy in the presence of a piagnetic fleld, magnetotun- 
neling, is a promising methodcj. In the present paper, we 
apply our idea to models of the above mentioned quasi- 
2D superconductors, namely the high-Tc cuprates and 
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K-{ET)2X. 

As regards k-(ET)2X, 
the pair potentiaU-n 
ous experimentstaEZH 



the existence of nodes of 
n suggested from 

However, the posi 
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tion of nodes is still a controversial issi. 
early theories . support d^2_y2-wave paLung 
experimentaHaa support dxy pairing^J. Stimulated by 
these recent experiments, two of the present authors re- 
visited this issue theoretically, and found that a d^y-^i^e 
pairing can slightly dominate over dj.2_y2 pairing when 
the dipierization of the BEDT-TTF molecules is not so 
strongl3. Thus, it is intriguing to clarify what is ex- 
pected in the magnetotunneling spectroscopy for each of 
the plausible pairing symmetries. 

We calculate surface density of states (SDOS) at the 
Fermi energy, which corresponds to zeijo-bias tunneling 
conductance in the high barrier limitu. We find that 
SDOS oscillates as a function of the direction of the mag- 
netic field, and has a broad maximum for the magnetic 
field applied in the antinodal direction, and a minimum 
for the magnetic field along the node. Although this re- 
sult is similar to those obtained in Ref. ^ where the bulk 
density of states is considered for free electrons (namely 
for a round Fermi surface) , it is by no means evident from 
the beginning whether similar results can be obtained at 
surfaces and for realistic lattice structures and band fill- 
ings (namely for realistic shapes of the Fermi surface). 
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FIG. 1: (a) Magnetic field H rotated with the angle d in xy 
plane, (b) Square lattice and (c) anisotropic triangular lattice 
with next- nearest neighbor hopping i! . Solid lines denote the 
Cooper pairs in real space. 



A schematic illustration of the models of quasi-2D su- 
perconductors and the direction of the magnetic field are 
shown in Fig. ^ where magnetic field is rotated within a 
xy-plane. Since the magnitude of the penetration depth 
is much more larger than that of coherence length in 
both superconductors, the vector potential can be ex- 
pressed to be A(r) = {HXe''^^ sind, -HXe^^^ cos6,0), 



where 6 is the direction of the magnetic field as defined in 
Fig. |l|. Thus the quasiparticle momenta kx and ky in the 



van- X and y directions can be chosen as kx = kx 



H 



COS 6 

and ky = ky — -j^^j-sinO, where He — <po / {'^'^ with 
<^o = V(2e). 

Within the BCS mean-field theory, in terms of the 
eigenenergy E and the wave functions Wfc, Vk, the 
Bogoliubov-de Gennes equation with the magnetic field 
applied parallel to the xy-plane is given by 



Vk 



= E 



Uk 
Vk 



(1) 
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{ik - £,-k) ± V&+?-fc)' + 4|Afc|2 , (2) 
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(3) 



with rik=£,k+ C-fc and Ffe = ^ttJ+IA^. 

In the following, we discuss the k dependence of ^k 
and Afc for the two cases considered, (i) the cuprates: 
We adopt the extended t-J model on 2D square lattice, 
which is considered to describe the low-energy excita- 
tions as a function of doping concentration S and superex- 
change interaction J. The reason why we use this model 
is that we can systematically investigate the doping de- 
pendence of Fermi surface. Since it is difficult to treat 
the constraint in the t-J model, we apply the Gutzwiller 
approximatioiO. In this case, and in Eq.(|l|) are 
given by 

^k = -2 (tcff + Jeff x) (cos kxa + cos kyoj 

—At'^g cos kxa cos kya — (4) 

and 



Afc = 2Ao (cos kxU — cos kyo) . 



(5) 



with tcff — Y+ji, and Joff = | (i+sy^ ^ ■ Here, the quan- 
tities X, /i, and Ag are determined self-consistently for 
each doping ratio 5. In the following calculation, we fix 
J /t — 0.3 and t' /t — —0.5 \t' /t — 0.5] as a typical values 
of hole-doped [electron-doped] cases. 

(ii) K-(ET)2A': Although this material has four 
molecules per unit cell (four bands), we mayjieduce it 
to a single band model for the present purposecJCJ. The 
single band description is not sufficient foi-determining 
which pairing symmetry is more favorable J13 but it suf- 
fices for the argument here, where we discuss the tunnel- 
ing spectrum for a given pairing symmetry. In the single 
band description, the quantity in Eq.(pl) is given by 



^fc — '-2t [coskxCL + coskya 



-2t' cos{kx + ky)a — fi. 



(6) 
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The values of t, t', and /z are chosen as to repro- 
duce the FecHii surface observed by Shubnikov-de Haas 
experimentaiS. As for plausible pairing symmetries in n- 
(EThX, we consider the d^2_y2-weLve pairing given by 
Eq.(||), and a d^y-like pairing given by 

Afc = 2Ao [cos k^a + cos kyU — a cos{k^a + kyu)] , (7) 

with a ~ 0.8. These pairing symmetries have been found 
to closely compete with each other in Ref. 

In order to compare our theory with scanning tunnel- 
ing microscopy (STM) experiments, we assume that the 
STM tip is metallic with a constant density of states, 
and that the tunneling occurs only for the site nearest 
to the tip. This has been shown to be valid through 
the study of tunneling conductance of unconventional 
superconductorso. The tunneling conductance spectrum 
at zero-energy is then given at low temperatures by the 
normalized SD093, 



t'/t = -0.5 



t'/t = 0.5 
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2kBT, 

k 

+ \vkf [S{uj + Ek) + d{uj + E^k)]} ■ 



(9) 



Here ps{(^) denotes the SDOS for the superconducting 
state while pn(^) the bulk density of states in the normal 
state. 
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FIG. 2: The angle variation of SDOS for the model of 
cuprates in ii" = 0.5Hc. (a) t'/t = -0.5 and (b) 0.5. 



First, we look into the case of the cuprates. p{9) is 
plotted in Fig. |^ for various 5. In the overdoped regime, 
SDOS exhibits clear fourfold oscillations, which has a 
maximum for the magnetic field applied in the antinodal 
direction, and a minimum for the magnetic field along 
the nodes. Thus, through the analysis of the position of 
maxima and minima of p{9), we can identify the nodal 
and antinodal directions. With the increase of (5, the 
magnitude of the difference between the maximum and 
minimum is enhanced. In the underdoped region, the 








FIG. 3: The Fermi surface and d-wave pairings. Dots lines 
represent the nodal lines, (a) a model of cuprates with a 
dj;2_j^2-wave: left and right panels are plotted for t'/t — —0.5 
and 0.5, respectively, (b) a model of k-(ET)2X for d^2_y2- 
wave (left) and da;j,-like (right). 



difference of the magnitude between maximum and min- 
imum is reduced reficcting on the square shape of Fermi 
surface (see Fig. |^). We now move on to the case of 
K-(ET)2X . For both of d^2_y2 and dxy-\i)^e cases, p{9) 
exhibits a characteristic oscillation reflecting the nodal 
structure and the twofold symmetry of the Fermi surface. 
For d3.2_j,2-wave pairing, p{6) has local minimum [max- 
imum] around (2n -I- 1)71 /A [wk /2] with integer n, while 
for dxy-wave pairing, p{0) has local minimum [maximum] 
around n7r/2 [(2n -|- l)7r/4]. As seen from this result, we 
can discriminate da;2_y2-wave from d^y 

To summarize, we have proposed tunneling spec- 
troscopy in the presence of a magnetic field, magnetotun- 
neling, for quasi-2D anisotropic superconductors, where 
the magnetic field is rotated within the 2D plane. SDOS 
at the Fermi energy, p{9), exhibits a characteristic oscil- 
lation upon rotating the direction of the magnetic field. 
As case studies, we have applied this idea to the high Tc 
cuprates and k-(ET)2X salts. p{9) is found to take its 
(local) minimum value when the applied magnetic field is 
parallel to the nodal direction of the pair potential inde- 
pendent of the detailed shape of the Fermi surface. The 
origin of the oscillation p{9) is considered to be essen- 
tially the same as that given in Ref. ^ but we stress 
that it is not expected from the beginning that the phase 
of the oscillation is the same as that in Ref. ^ for ar- 
bitrary doping concentrations and/or for various lattice 
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FIG. 4: The angle variation of SDOS for a model of k- 
(ET)2X. 



structures. In fact, we have seen that in the under- 
doped regime of the cuprates, the oscillation, although it 
still has the same phase as that in Ref. becomes ex- 
tremely weak reflecting the characteristic band structure 
near half-filling. Although only two kinds of quasi-2D 
superconductors are studied here, the proposed magne- 
totunneling spectroscopy should serve as a strong probe 
to identify the pairing symmetry for various unconven- 
tional superconductors. 
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